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Molecular spiders are synthetic molecular motors based on DNA nanotechnology. While natural 
molecular motors have evolved towards very high efficiency, it remains a major challenge to develop 
efficient designs for man-made molecular motors. Inspired by biological motor proteins like kinesin 
and myosin, molecular spiders comprise a body and several legs. The legs walk on a lattice that is 
coated with substrate which can be cleaved catalytically. We propose a novel molecular spider design 
in which n spiders form a team. Our theoretical considerations show that coupling several spiders 
together alters the dynamics of the resulting team significantly. Although spiders operate at a scale 
where diffusion is dominant, spider teams can be tuned to behave nearly ballistic, which results 
in fast and predictable motion. Based on the separation of time scales of substrate and product 
dwell times, we develop a theory which utilises equivalence classes to coarse-grain the micro-state 
space. In addition, we calculate diffusion coefficients of the spider teams, employing a mapping 
of an n-spider team to an n-dimensional random walker on a confined lattice. We validate these 
results with Monte Carlo simulations and predict optimal parameters of the molecular spider team 
architecture which makes their motion most directed and max;imally predictable. 

PACS numbers: 87.16.Nn, 82.39.Fk, 05.40.Fb, 02.50.Ey 



I. INTRODUCTION 

How the motion of molecules along predefined traf- 
fic routes emerges and how these molecules self-organise 
is now an experimentally tractable question due to ad- 
vances in nanotechnology. Molecular motors that have 
evolved inside cells and perform well-defined tasks [1] in- 
spired the engineering of DNA devices performing motor 
business on the nanoscale [2-4]: so-called DNA walkers 
have been build that move or diffuse along a substrate [5- 
7]. Among the first autonomous synthetic walkers was 
a motor design that used a catalytic reaction to cleave 
a substrate in order to move forward [8]. Since then, 
a plethora of different motor molecules have been built 
from scratch in the laboratory. They do not only serve 
technological advances, but also shed light on the ba- 
sic principles of molecular movement, e.g. of biological 
molecular motors. One class of molecules that attracted 
a great deal of attention are molecular spiders [9] . They 
combine the catalytic activity of nucleic acids with a mul- 
tivalent design: attached to a body are several legs of 
single-stranded DNA. These DNA legs can bind to and 
catalytically cleave a substrate. This can be repeated 
over and over again, which in turn generates processive 
motion: while individual legs dissociate from the sub- 
strate on a timescale of seconds, the multipedal architec- 
ture ensures tight binding of the spider to the substrate 
for hours [9]. Recent experiments used DNA origami 
to build quasi-one-dimensional tracks for molecular spi- 
ders [10]. A predescribed substrate landscape allows to 
assign special tasks to a spider and for instance control 
its movement. The simple yet well-defined design makes 
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it possible to study spiders in great detail and probe the- 
oretical predictions. 

Molecular spiders have also been theoretically studied 
extensively in recent years. Antal et al. [11] and Antal 
and Krapivsky [12] were the first to propose an abstract 
model that describes the dynamics of molecular spiders. 
They analysed the spiders' kinetics for various architec- 
tures and found a variety of interesting effects which arise 
due to the mutual exclusion of spider legs on the lat- 
tice and the presence of the substrate. Substrates are 
cleaved slowly in comparison to hopping from already 
cleaved sites. This distinction leads to subtle memory 
effects that affect the spiders dynamics and result in a 
bias towards the substrate [12]. When the spider is in an 
all-cleaved area, principles emerging from simple exclu- 
sion processes [13, 14] allow a derivation of the spiders' 
diffusion constants [15]. 

In the meantime, mechanistically more detailed sys- 
tems have been considered. These include the varia- 
tion of the rate constants involved in the chemical reac- 
tions [16, 17] and boundary conditions [17], as well as the 
number and length of legs [18]. Samii et al. [17] investi- 
gated the spiders' stepping gait and considered inchworm 
as well as hand-over-hand spiders. Semenov et al. [16] 
showed that spiders experience a rather extended time 
period of superdiffusion given that the cleavage rate r is 
small. More complex spiders in quasi-one [19] and in two 
dimensions [20] have also been studied. Moreover, there 
have also been recent studies focussing on mathematical 
aspects like recurrence, transience and ergodicity [21, 22], 
as well as random environments [23, 24]. These investi- 
gations have examined molecular spiders independently 
from their chemical motivation as a general class of mul- 
tivalent random walkers [19]. 

The rich variety and diversity of these recent stud- 
ies show that molecular spiders are a versatile system 
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to study artificial molecular motors both theoretically as 
well as experimentally. However, many challenges still 
remain in improving their efficiency and tailoring the spi- 
ders' design for possible biotechnological applications [7] . 

In this study, we examine dynamic and stochastic 
properties of a novel molecular spider design: n molecu- 
lar spiders are constrained due to their joint attachment 
to a single linking node which may be considered as a 
primitive model of a cargo. The resulting spider-spider 
interactions lead to collective effects which enhance the 
motor properties of the n-spider team. We show that 
spider teams are faster and move more persistently along 
their track than individual spiders. We also predict that 
the spider teams move at reduced randomness and thus 
are candidates for applications that require reliable, i.e. 
predictable motion [4]. 

This paper is organised as follows: in Sec. II we pro- 
vide a detailed picture of how molecular spiders function 
and give a comprehensive introduction to the existing 
theoretical models before we define the dynamics of an 
n-spider team. Subsequently, in Sec. Ill A we present our 
main results: spider teams have enhanced motor proper- 
ties. To explain these numerical observations, we give 
a comprehensive analysis of the stochastic dynamics of 
a spider team. In particular, we perform a reduction of 
the state space of the spider teams and thereby calculate 
the mean number of consecutive directed steps a spider 
team performs while attached to the substrate bound- 
ary (Sec. IIIB). Moreover, we explore the validity of the 
resulting network representation of the spider team dy- 
namics and also show how it breaks down (Sec. Ill C). In 
addition to this approach, we provide an exact mapping 
of the n-spider team to an n-dimcnsional confined ran- 
dom walk (Sec. HID). This enables us to quantify the 
diffusion coefficient which describes the motion of a spi- 
der team during diffusive periods (Sec. HIE). Finally, in 
Sec. IV we bridge theoretical and experimental observ- 
ables and predict the existence of optimal parameters 
which maximise the spider teams' predictability. Finally 
(Sec. V), we conclude and identify connections to related 
fields. 



II. MODEL DEFINITION 

Our model is based on the theoretical description of 
molecular spiders introduced by Antal et al. [11] and 
Antal and Krapivsky [12] that was motivated by exper- 
iments of Pei et al. [9]. They propose a spider design 
that consists of a central body and I legs that are at- 
tached to it. Each leg has a certain length and thus the 
overall spider can span a maximal distance s. In the 
experiment, a spider is exposed to a (one-dimensional) 
lattice, to which a substrate is attached. Since binding 
of leg and substrate happens through the Watson-Crick 
mechanism [25], only one leg may bind to a lattice site 
at a time. In the model, this corresponds to an exclu- 
sion process in that the movement of one spider leg is 
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FIG. 1. (Colour online) Cartoon of the spider team model 
and definition of the leash length d. (a) Two spiders are at- 
tached to a joint cargo with an inelastic string. Both spiders 
walk on their respective one-dimensional track. Hats indicate 
the presence of substrate, (b) The finite length of the link- 
ing string induces a maximal distance between the spiders' 
bodies which gives rise to a maximal span of the spider team, 
characterised by the "leash length" d. 



constrained by the spider's remaining legs. The lattice 
prevails in two states: with and without substrate. Legs 
which bind to lattice sites with substrate can remove it 
(chemically: they cleave it, only a shorter part remains 
bound to the lattice) , which happens along with unbind- 
ing from that site at rate r. By contrast, spiders unbind 
from sites without substrate {i.e. from product sites) 
at rate 1. In the model, a substrate is always cleaved 
when a leg steps away from it, and rebinding of a leg 
to a new lattice site happens instantaneously. Two dif- 
ferent rules to rebind to a new lattice site have to be 
distinguished: spiders' legs either have a certain order- 
ing, i.e. they cannot "overtake" each other; these spiders 
are termed inchworm spiders [11, 12, 16]. Alternatively, 
spider legs have no ordering, they can step over each 
other; those spiders have been called "quick spiders" [11] 
or "hand-over-hand" spiders [17, 18] in previous studies. 
Both types of spiders show quite different behaviour [18] 
and have to be well distinguished. In this paper, we will 
concentrate on inchworm spiders. 

Although in our model a leg which has just unbound 
from the lattice rebinds to the lattice instantaneously, we 
allow a spider's leg to rebind to any lattice site as long as 
the new leg configuration does not violate any of the re- 
strictions imposed by the leg length or the ordering of the 
legs. In particular, this implies that rebinding to the lat- 
tice site from which the leg just unbound is possible [26]; 
this can be motivated from experiments where the typi- 
cal timescales for binding to substrates exceed those for 
diffusion by orders of magnitude [17]. In addition, our 
choice obviates unphysical situations that might occur 
for spider teams due to the complete blockage of a leg. 

Hollow circles (o) denote unoccupied lattice sites, filled 
circles (•) indicate that a leg is attached to that site. 
The presence of substrate is marked with a hat, i.e. •, 
or o. Throughout this paper, we consider bipedal spiders 
(i.e. I — 2) with a maximal leg span of s = 2. Spiders 
may thus only arise in either the spanned (•o») or the 
relaxed (••) configuration. For this case, the geometry 
of the cleaved sites, which is usually called product sea, 
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is an interval on the one-dimensional lattice; it gives rise 
to memory effects which stem from irreversible substrate 
cleavage [11]. 

Samii et al. [18] suggested that the lattice could be 
prepared with substrates on the right, and products on 
the left hand side from the very beginning, and called this 
initial condition P-S lattice. This asymmetry makes some 
calculations easier, and it provides a symmetry breaking 
direction already at the beginning of the dynamics. We 
are going to use this kind of lattice throughout this paper. 

Taken together, the spiders which we examine in this 
study are bipedal {I = 2) inchworm spiders with a maxi- 
mal span of s = 2, which walk on a one-dimensional P-S 
lattice. Every spider's leg may rebind to any accessible 
lattice site as long as the ordering is preserved, including 
the site from where it just unbound. 

Based on this model for molecular spiders, we pro- 
pose a minimal model for a team of molecular spiders. 
Several, say n, molecular spiders are linked to a (vir- 
tual) cargo with an inelastic leash (i.e. a string) of a 
well-defined length d. Each of these spiders runs on its 
own one-dimensional track. This is similar to biological 
molecular motors like kinesin-1 [27, 28] that walk along 
one-dimensional microtubule filaments [1] . We call these 
ensembles of spiders that jointly pull a cargo a spider 
team. For a cartoon of a team of two spiders, see Fig. 1. 

Note that the role of the "cargo" is not primarily to 
put load on the spiders, actually we set the mass of the 
cargo equal to zero. In contrast, the cargo mediates the 
interaction among the n spiders comprising the team: 
since the strings used for linking the spiders to the cargo 
are inelastic with some length a, any two of the spiders' 
bodies may mostly be 2a away from each other. From the 
bodies, the furthestmost reachable lattice site is given by 
the spiders' legs' length, call it b, so that the maximal 
distance between the leftmost and the rightmost leg of 
all the spiders in the team is given by 2(a -I- b) =: d. 
Mathematically, letting [pi] denote the position of the 
ith. spider's left [right] leg, this restriction reads 

1p,-a,i <dVz,j. (1) 

Note that this is a global constraint which restricts the 
spider team, in contrast to the local constraint limiting 
the span of an individual spider, 

\p,-X,\ <s EE2Vi. (2) 

The definition of d is visualised for a 2-spider team in 
Fig. 1(b). For simplicity of language, and to capture an 
intuitive understanding especially for 2-spider teams, we 
will call d the leash length in the following. 

III. RESULTS 
A. Enhanced properties of n-spider teams 

We performed extensive numerical simulations to char- 
acterise the dynamic properties of n-spider teams. Our 



simulation data show that the constraint arising through 
the leash that holds the spider team together induces col- 
lective effects among the n spiders. We find that the in- 
corporation of a spider into a team enhances many of the 
motor properties: the mean travelled distance of a spider 
team exceeds that of single spiders by far, up to orders of 
magnitude, for a rather small cleavage rate r = 0.01, see 
Figs. 2(a), (b). In addition, a spider team's movement is 
a lot more "predictable" . This can be inferred from the 
width of the probability distributions, see Fig. 2(a), and 
the shaded areas depicted in Fig. 2(b), which illustrate 
the standard deviation of the mean displacement. 

Another important quantity is the mean square dis- 
placement (MSD) of the spider teams, see Fig. 2(c). It 
shows a steep increase at intermediate timescales, sim- 
ilar but stronger and longer-lasting compared to recent 
results by Semenov et al. [16] for single spiders: in this 
regime spiders move superdiffusively. To quantify the 
time-dependent effects of supcrdiffusion, we evaluated 
the "slope" of the variance in a double logarithmic scal- 
ing, i.e. the effective exponent 

rflog((x(t)-(x(t)))^) 

= = dl^t ' 

which provides a measure for diffusivity, see also 
Refs. [16, 19, 29]. Figure 2(d) shows a{t) for a single 
spider and several different spider teams. Remarkably, 
the 4-spider team travels almost ballistically and the pe- 
riods of "instantaneous supcrdiffusion" of spider teams 
(i.e. times with a > 1.1 [16]) last much longer compared 
to single spiders. The nontrivial shape of a{t) indicates 
the multitude of dynamic processes that are involved in 
the spider team's dynamics: initially, a ~ I for t < 1 
for all configurations, reflecting the very first hop of the 
spiders' left legs. In succession, until t < = 100, 
the spiders' right legs have typically not yet cleaved a 
substrate, whereas the left legs jump back and forth, 
hence the variance is approximately constant and thus 
a < 1 (for these two regimes, see also a more explicit 
discussion in Ref. [16]). Had we chosen other starting 
conditions for the spiders, the behaviour at short time 
would look different. Likewise, also the following regime 
until t < 10^ . . . lO'^ results from the fixed starting con- 
ditions: while at early times the spider team does not 
feel the leash and all spiders can move independent from 
each other, at some point the leash is fully spanned and 
the spiders at the most extreme position (i.e. those con- 
tributing most to the variance) are retarded. This leads 
to a transient decrease of a. This regime is unique to 
spider teams since it is an effect constituted by the leash. 
Finally, for large times t > 10^ . . . 10'^, the memory of ini- 
tial conditions is lost and a becomes maximal. Clearly, 
the maximal value of a is greatest for n — A of the dis- 
played configuration. As time increases further, a de- 
creases slowly which is due to the fact that more and more 
spiders move away from the product-substrate boundary 
(see also Ref. [16]). Figure 2(e) shows the velocity of 
the spider team by means of the derivative of the mean 
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FIG. 2. (Colour online) Dynamic properties of spider teams. Positions are given in lattice units throughout this work; time is 
defined by setting the hopping rate fi-om products to 1. Thin shaded lines show data from finite difference approximations, thick 
lines show smoothing Bezier curves, (a) Probability distributions (histograms) of spiders to be at position x at time t = 10^; 
simulation data were binned with a box size 1. Depicted are distributions for a single spider and spider teams with n = 2, 3, 4 
and d = 8, and cleavage rate r = 0.01. While the single spider distribution follows nearly a Gaussian centred close to the 
origin, the distributions of spider teams are clearly skewed and shifted towards larger x. The asymmetry stems from the P-S 
preparation of the lattice at t = (products at the left, substrates at the right) [17]. (b) Mean displacement as a function 
of time (lines). The shaded areas represent the standard deviation around the mean displacement for a single spider and the 
4-spider team, respectively, and provide a measure for the randomness of the spiders' motion. Note that the visual impression 
of the standard deviation is rather that of a relative deviation, since the plot is in double logarithmic scale. For a discussion 
of the temporal behaviour of the dynamic quantities see the main text, (c) Mean square displacement (MSD) as a function of 
time, {x^{t)). (d) The variance's effective exponent a[t), see Eq. (3). For diffusion, {{x{t) — {a;(t)))^) oc , hence a — 1; 
superdiffusion corresponds to a > 1.1 [16], and ballistic motion to a = 2. The superdiffusive regime of spider teams lasts longer 
than that of single spiders; large spider teams reach nearly-ballistic motion for significantly long times, (e) Mean velocity of the 
spiders as a function of time. The mean velocity is defined as the time-derivative of the mean displacement, d{x{t)) /dt. Spider 
teams outperform single spiders by an order of magnitude, (f) Sample trajectory of a single spider (top), and a 4-spider team 
with d = 8 (bottom). Periods in which the spider (team) is in the vicinity of the product-substrate boundary are shaded. 



displacement with respect to time. Clearly, the velocity 
of a 4-spider team outperforms that of a single spider by 
more than one order of magnitude. 

These pronounced effects are in a way surprising: at 
first sight, one might speculate that the coupling leash 
which imposes an additional constraint on the spiders 
would handicap the spider team's motion and make it 
slower. This is clearly not the case. To the contrary, the 
dynamic properties of the spider teams are enhanced. In 
the remainder of this section we will explain this effect 
using analytical arguments. 



B. Boundary periods 

a. Single spiders. Key to the understanding of an 
individual molecular spider's motion is to unravel the 
mechanism for biased motion. To this end we distin- 
guish between two qualitatively different dynamic states 
of the spiders: looking at single trajectories of molecular 
spiders we find that there are periods of time in which 
the spider's motion is strongly directed, and other peri- 
ods with undirected, diffusive motion (see Fig. 2(f)). In 
the following, we will call these dynamic states boundary 
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FIG. 3. (Colour online) Definition of a boundary period, (a) Path of a single spider through a boundary period. The period 
always starts in state (i). From there, the spider can change to (ii), and back. When the right leg cleaves the substrate, the 
spider arrives at {iv), (v), {vi), or {vii). Arriving at {vii) corresponds to continuing the same boundary period from a new 
substrate (with '^{vii) being the new (i)"), since {vii) and (i) are equivalent up to translation. Hence, the number of steps is 
raised by 1 upon arriving at {vii). If, by contrast, the spider reaches {iii), the boundary period ends and a diffusive period 
begins. The probability to make a successful step, i.e. to reach {vii) before {iii), is the bias p+ calculated by Antal and 
Krapivsky [12]. The number of steps during a boundary period is then the number of transitions {i) — >■ {vii), without reaching 
{iii) in between. This is equivalent to the number of cleavages during a boundary period, not counting the very last cleavage 
(which is not counted since by definition the spider steps away from the boundary after the last cleavage, and we only count 
forward steps), (b) Example of a boundary period of a 2-spider team, (a) None of the spiders is in a boundary period, hence 
none of them experiences a bias. Thus, the spider team is in a diffusive period. When the lower spider reaches a substrate (/3) 
it enters a boundary period. Thus, also the spider team enters a boundary period. In succession, the lower spider's right leg 
happens to cleave the substrate (7). The lower spider can then find its way to a new substrate {S) what constitutes | successful 
step for the spider team and preserves the boundary period. If the upper spider, in this case, steps to a substrate (e), this does, 
however, not yet constitute a step. This is because although the spider team is in a boundary period, the upper spider has not 
been in a boundary period itself during this team's boundary period. Since a step essentially reflects a cleavage, no step can be 
integrated in this case. If the lower spider steps away from the new substrate (C), the spider team enters a diffusive period. In 
analogy to single spiders, the number of steps during a spider team's boundary period is equivalent to the number of cleavages 
during that period, divided by the number of spiders, and not counting each spider's last cleavage event. 



periods and dijfusive periods, respectively. To define the 
notion of these periods, it is convenient to distinguish 
between the steps of the spider's legs and the step of the 
spider as a whole. We define a spider step as a tran- 
sition from a spread configuration (o»o»o) to another 
spread configuration shifted by one lattice unit forwards 
or backwards, i.e. oo»o»or»o»oo, irrespective of the 
sites being products or substrates. During a diffusive pe- 
riod all the spider's legs are attached to product sites and 
therefore the spider steps with equal probability in both 
directions [11]. In contrast, biased spider motion can 
emerge in the vicinity of the boundary between prod- 
uct and substrate sites. We define a boundary period as 
follows: it starts with a spread configuration where the 
right spider leg is attached to a substrate (•••o»oio--.), 
and ends when the spider has fully stepped away from 
the substrate boundary (••••o»o5---), (the dots indi- 
cate that the block of displayed lattice sites may have 
been shifted during the boundary period) as illustrated 
in Fig. 3(a). As a consequence, during a boundary period 
the substrate boundary is shifted by an integer number 
of lattice units forward. 

For single spiders the bias can be measured by cal- 
culating the first passage probability, p+, for the spider 
to progress one step forward during a boundary period, 
i.e. p+ =Prob{»o»oo o»o»o ^ oo»o4}; see also 
Fig. 3(a) for an illustration of the corresponding dy- 
namic processes. By analysing all possible sequences of 
transitions, Antal and Krapivsky found an explicit ex- 
pression for the bias, namely p+{r) = [12], valid 



for spiders with legs always jumping to neighbouring 
sites. Similar calculations can be performed for spi- 
ders whose legs may also rebind to the same site again 
(like those we consider throughout this paper), leading 
to p+{r) — g^g^ . The mathematical expressions for 
p+ and p+ differ only slightly, in particular they are 
equal in the limits lim^-i-o P+ ('') = linir— +oP+(^) — | E^nd 
p+(r = l)=p+(r = l) = i [30]. 

There is a special feature of single spiders which makes 
the definition of pj^. straightforward in this case: the 
spread configuration •o» of the spider's legs is unique, 
since a spider step to the right correspond to a transla- 
tion of both legs to the right, and hence the configuration 
before and after a step is the same; cf. 3(a) (i), {iii), and 
{vii). As we will show below, this is a property which 
unfortunately does not extend to spider teams. 

A quantity which does not require this uniqueness is 
the mean number of consecutive directed steps that a spi- 
der performs during one boundary period. This quantity 
will be denoted (S) in the following. For single spiders, 
{S) can be calculated as 

00 

(5(p+))=^,p, = ^, (4) 

where 

p, = {p+)\l-p+) (5) 

is the probability that the spider walks precisely j steps 
during a boundary period, before it leaves the boundary 
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and enters a diffusive period. Let us emphasise that {S) 
is different from the mean "number of steps the spider 
makes in the B state" [16], (Sb), as defined by Semenov 
et al., which counts the number of leg movements ("leg 
steps" in our terminology). By contrast, (S) only counts 
a step if both legs have been shifted to the right without 
having moved to the left ( "spider steps" ) , i.e. the number 
of times the spider consecutively reaches (vii) before (Hi), 
starting from (i) in Fig. 3(a). 

The number of consecutive spider steps, (S), is equiva- 
lent to the number of cleavage events during a boundary 
period. Not counted is the last cleavage before the spider 
leaves the boundary period, since this corresponds to a 
backward step of the spider, cf. Eq. (5). 

b. Spider teams. Clearly, the motion of a single spi- 
der is biased only during boundary periods, and undi- 
rected during diffusive periods. However, it is manifest 
that a spider team's motion is not completely diffusive 
as long as any of the spiders comprising the team is in 
a boundary period. Hence, we consider the spider team 
being in a boundary period if at least one of its spiders re- 
sides in a boundary period. In order to compare the per- 
formance of individual spiders with that of spider teams, 
it is now essential to find a way how to count the number 
of a spider team's steps during a boundary period. Basi- 
cally, a team moves forward by one step if the boundary 
between substrate and product sites is shifted forward 
by one lattice unit on average. To this end we count 
every cleavage event but for each spider's last cleavage 
before the team leaves the boundary period. In analogy 
to a single spider, the latter avoids counting those events 
where the spider team moves away from the boundary 
and thereby steps backward, cf. Fig. 3(b). The number 
of steps of a spider team is then given by the number of 
such cleavage events divided by the number of spiders in 
a team, in analogy to fractional steps of molecular motors 
like kinesin [31]. For example. 



o»o»o oo»»o 

„ . . — > 

o»»oo oo»o» 



(6) 



corresponds to two steps of the lower spider and thus one 
step for the spider team. 

As we consider two or more coupled spiders, the trans- 
lational symmetry of the state before and after a com- 
plete step (•o46 and o»o», respectively for a single spi- 
der) is broken, likewise the uniqueness of the state which 
is the first during a boundary period (•oi for a single 
spider), is lost. For example, 
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all are possible states at the beginning of a boundary pe- 
riod. It is therefore no longer possible to calculate the 
probability to step to the right (denoted p+ for single spi- 
ders) without further specification of these initial states. 
For spider teams the probability for a forward step ex- 
plicitly depends on the particular state from which it 
starts. 



This complexity prohibits an analytic treatment of the 
stochastic dynamics in general. However, if the relative 
rate of substrate cleavage is small compared to the rate 
of hopping from product sites, r ^ 1, the dynamics be- 
come amenable to a theoretical analysis. While in this 
limit the motion of the boundary between substrate and 
product sites is slow, the dynamics of spider legs bound 
to product sites are fast. This suggests to group states 
into classes characterised by the slow variable, i.e. the 
distance between the ends of the product seas, denoted 
by A. In addition, it turns out to be convenient to in- 
troduce subclasses according to the number of spiders 
attached to substrates, a. In the following we will illus- 
trate this for teams comprised of rt = 2 spiders and a 
leash length d = 2. All states 
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comprise the class 
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(8) 
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Likewise, configurations with A = and a = 1, i.e. with 
only one spider having a leg at the boundary, are possible: 



o»o»o 
o»»oo 



- [Oi]. 



(10) 



Here, we made use of the invariance under renumbering 
of spiders, it is irrelevant if we label the "upper" spider as 
1 and the "lower" as 2, or the other way round. Hence, 
irrespective of whether the lower or the upper spider's 
leg is bound to a substrate, both contribute to class [Oi] . 
That same renumbering symmetry can also be applied 
when one considers states where the lower and the upper 
product seas do not end at the same position. This leads 
to the classes 



o • o • o 
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=: [I2] and 
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as well as 
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[2i 



=:[li], (11) 



(12) 



This completes the list of possible classes with cr 7^ since 
the constraint on the leash length forbids class [22] , as 
well as classes [A^.] with A > 2. For general d, class [c?2] 
and classes with A > d are not allowed. 

One can show that the classification of states by means 
of the distance of the product seas' ends and the num- 
ber of spiders at the boundary is reflexive, symmetric, 
and transitive, and hence defines an equivalence relation. 
Therefore, we tentatively used the symbols ~ and [ • ] 
in the previous equations. 

Instead of a large number of "micro" -states, we are 
now left with only five equivalence classes which include 
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all the spider states at the boundary. The reduction of 
complexity can be pushed even further: classes [Ai] with 
only one leg attached to the substrate are only transient 
in the sense that they will always decay into classes with 
two legs attached [A2] (as long as A < rf). Consider, for 
example, a spider team in class [Oi] where one spider's 
right leg is attached to a substrate while the other spi- 
der's legs are free to move on product sites. Since the 
diffusion time of legs on products is small compared to 
the expected residence time l/r of the leg on the sub- 
strate, the transition [Oi] — > [O2] is almost certain and 
happens on a time scale ~ 1 (fast compared to substrate 
cleavage). 

All possible transitions between the classes can be vi- 
sualised as the following reaction scheme: 



[Oi] [li] 
[O2] ^[12] ^■ 



[id -1)0] 
i 

[(d-l)i] 



diffusive 
period 



(13) 



where A is constant along a column and a along a row, re- 
spectively. As explained above, vertical transitions from 
[Ai] to [A2] are fast (emphasised with bold red/grey 
arrows in Eq. (13)). In contrast, horizontal and diagonal 
transitions involving substrate cleaving events and hence 
leading to A — )■ A ± 1 are slow. Since vertical transitions 
occur with certainty and fast, we can eliminate the tran- 
sient classes [Ai] and reduce to a reaction scheme for 
the most stable subclass of each class, shown in boldface 
in Eq. (13) and signified [A] in the following: 



[0] 4^ [1] 



-n diffusive , , 
period ' ^ ' 



The numbers above and below the arrows are transition 
probabilities into the respective classes, reflecting that 
each of the two spiders may cleave a substrate with equal 
probability for A < d. The class [d] has to be treated 
separately as it constitutes a gate from the boundary into 
the diffusive period. 

Our next set of tasks is now threefold: first, in order 
for our classification scheme to be a consistent reduction 
of the stochastic processes, all states comprising the gate 
class [d] = [di] should have the same survival probabil- 
ity n, i.e. the same probability not to exit into a diffu- 
sive period. This is indeed the case for sufficiently small 
cleavage rates r: in the limit r 0, substrate cleavage 
events are rare compared to hopping from product sites. 
Therefore, the dynamics exhibit a time scale separation 
where all the legs attached to products quickly visit any 
accessible lattice site while the legs on substrate sites re- 
main stuck. In other words, the dynamics within class 
[di] are ergodic and equilibrate, and all "micro" -states 
effectively reduce to one coarse-grained "macro"-state, 
namely the class [c?i] . Second, we have to calculate the 
survival probability n by analysing all the various routes 



between the "micro" -states. Third, in order to determine 
the mean number of consecutive steps (S), the reduced 
reaction scheme of Eq. (14) has to be solved. 

We now address the calculation of the survival prob- 
ability n. In principle, this can be done for arbitrary 
complex spider teams. For the purpose of illustration, 
we continue the example from above with 2 spiders and 
a leash length d = 2. We consider all states comprising 
class \2i \ . These are 
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and their respective "mirrored" states, i.e. the states 
with spider 1 and 2 interchanged. Let us illustrate the 
calculation for the particular initial state 0. Legs un- 
bind from products at rate 1 and from substrates with 
rate r. Hence, from this configuration, the probability 
that the upper right, or the lower right leg is the first 
one to unbind is 1/(3 -I- r), and r/(3 -I- r), respectively. 
The left legs unbind first with probability 1/(3 -I- r) each. 
If now, for instance, the lower right leg detaches, it may 
either reattach to the very same lattice site again, or it 
may step one site to the right. In either case it cleaves a 
substrate. Both processes happen with equal probability. 
Hence, altogether, the transition probability for the lower 
right leg to step to the right is given by 7'/2(3 -I- r). The 
analysis can be continued from the resulting states, and 
finally leads to a high dimensional system of linear equa- 
tions. The results obtained by solving the ensuing sets 
of equations are shown in Fig. 4 for all initial states com- 
prising class [c?i] . Clearly, as r approaches 0, all survival 
probabilities approach a single value 

115 
176 

This result is reassuring, as it confirms our heuristic ar- 
guments on the equilibration of states within class [di] , 
and thereby justifies to combine several different states 
into one class in the limit r — )■ 0. 

All the complexity of calculating the mean number of 
steps {S) of a spider team during a boundary period has 
now been reduced to analysing the various routes be- 
tween the equivalence classes. Since each transition [32] 
in Eq. (14) corresponds to a directed step done during a 
boundary period, the number of these steps {S) is equiv- 
alent to the number of (undirected) jumps performed by 
a simple random walker with reflective, and absorbing 
boundary conditions on the left, and right end of the re- 
action scheme, respectively. As detailed in Appendix A, 
the general solution for the mean number of steps during 
a boundary period in the limit r — )■ 0, and for arbitrary 
d, reads 



H 



0.65. 



(16) 



{S[d,r-^ 0)) = ^ + (d - 1)^. (17) 
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FIG. 4. (Colour online) Justification for the equivalence 
classes in the limit r — >■ 0. Shown are the analytically cal- 
culated probabilities that a spider team (n = d = 2) success- 
fully completes one step during a boundary period, starting 
from the specific states ® - (B) as given in Eq. (15). Each 
line corresponds to a state of the equivalence class [2i] [cf. 
Eq. (12)]. In the limit r — ^ 0, the probability to step forward 
for all eight states collapses to a fixed value 11 ~ 0.65. 



FIG. 5. (Colour online) Validity of the equivalence class for- 
malism. Shown are the simulation results for the mean num- 
ber of steps, (S), for a 2-spider team and different values 
of d and cleavage rates r = 0.0001,0.001,0.01; broken lines 
are a guide to the eye. The theoretical result derived within 
the equivalence class formalism for r — >■ (black) is exact 
for d = 2,3,4 (Tab. I), and we assumed 11 = | for d > 5 
(Eq. (17)). 



For our example of a two-spider team with d 
obtain using Eq. (16) 
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0)3=1) 



291 
~61 



4.77. 



(18) 



We also analysed more complex spider teams with size 
n = 2, 3 and up to a leash length of d = 4, and found 
even larger mean step numbers, compared to | for a sin- 
gle spider. Obviously, during boundary periods even the 
simplest spider teams behave significantly more directed 
and progress a lot further on average, compared to indi- 
vidual spiders. This result is remarkable since directed 
motion is desirable for applications and a rare feature at 
the nanoscale. 

The analytical results are summarised in Tab. I where 
they are also compared with Monte Carlo simulations 
which match them at a very high accuracy. 



{S{r 0)2), analytic {S{r < lO"'')^), sim. 



n = 1 
n = 2,d = 2 
n = 2,d = 3 
n = 2,d = 4 
n = 3,d = 2 
n = 3, d = 3 



5 ^ 
3 ~ 
291 „ 
61 " 
3170931 , 
443341 ' 
4055316673 
414459263 
340881 „ 
48391 " 
16.3745... , 
1.34258... ' 



1.6667 
:i 4.770 
« 7.152 
Si 9.785 
7.044 
12.196 



1.6672 ±0.0015 
4.769 ±0.003 
7.146 ±0.005 
9.785 ± 0.008 
7.042 ± 0.006 
12.204 ±0.012 



TABLE I. Comparison of analytic and simulation results for 
the mean number of steps during a boundary period, (S). 
Analytic values were derived in the limit r — ^ 0, simulation 
results were obtained for very small r < lO"*. Simulations 
and analytical calculations show excellent agreement. 



C. Validity of the equivalence classes 



With increasing d, the spiders forming a team become 
more and more independent since it is increasingly un- 
likely that a spider "feels" the constraint of its team- 
mates. In particular, the probability 11 that a spider in 
class [di] reaches [(d— 1)2] without exiting the boundary 
period (cf. Eq. (14)), converges towards the probability 
p+ that a single spider makes a step to the right which is 
I for r -> 0. Hence, assuming H = | for large d, Eq. (17) 
would imply that the mean number of steps increases lin- 
early with d. Indeed, in the asymptotic limit r — >■ this 
agrees well with the simulation data. However, with in- 
creasing r deviations from this linear behaviour become 
more and more significant; cf. Fig. 5. 



This can be explained as follows: for increasing leash 
length d, the configuration space accessible to the spider 
team becomes progressively larger, so that it takes longer 
to completely exploit it, i.e. the equilibration time grows. 
Conversely, the average time of substrate cleavage scales 
as l/r. With increasing r and/or d these two timescales 
become comparable. The assumption of time scale sepa- 
ration, on which the reduction of the dynamics to equiv- 
alence classes was based on, then becomes invalid. Con- 
cluding, the equivalence class concept which we derived 
in the previous sections provides a very good approxi- 
mation for small but finite substrate cleavage rates r, as 
long as the leash length d is not too large. 
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FIG. 6. (Colour online) A spider team can be mapped to 
a random walk in a confined environment: transitions of a 
spider's leg correspond to a change of its centre of mass co- 
ordinate Ci of ±|. Shown is the mapping of a 2-spider team 
with a leash length d = 2. The shape of the environment 
(solid) follows from the leash constraint which confines the 
span of the spider team. From d = 2 follows that the left- 
most left and the rightmost right leg of the two spiders may 
be at most two lattice sites apart. With that restriction, the 
allowed configurations of the team follow directly, as can be 
seen with some explicit configurations in the left and the top 
part of the figure. 



D. An exact mapping to a confined random walker 

For a bipedal spider with a maximal span of s = 2, a 
single coordinate, the "centre of mass" -coordinate fully 
describes the position of the spider's legs. Hence, it is 
possible to map the motion of the single spider's legs on 
iZ, the set of integers and half-integers, with hopping of 
the legs corresponding to changes of the centre of mass. 
This mapping can be extended for a spider team: the po- 
sition of a n-spider team is characterised by a position on 
a n-dimensional square lattice where each of the n axes 
corresponds to the centre of mass of one of the spiders 
comprising the team. The dynamics of a spider team 
then correspond a trajectory on that lattice. However, 
due to the leash constraint, not all sites on this lattice are 
accessible to the spider team. To illustrate this, let us fo- 
cus on a two-spider team with leash length d = 2. Fixing 
the first spider's centre of mass ci, the other spider's cen- 
tre of mass C2 is restricted to be near Ci due to the leash 
constraint, cf. Fig. 6. We have to distinguish between 
two cases. Spider 1 is either in a spread or a relaxed 
configuration, e.g. Ci = or ci = i, respectively. If it is 
in the spread configuration ci — 0, then the other spider 
may be in one of three configurations: C2 £ {—^,0, ^}. 
For the relaxed configuration ci = i, there are five con- 
figurations possible for the second spider: —5,0, |, 1, and 



|. Geometrically, this leads to a staircase shape for the 
accessible set of states. For arbitrary d, the step width 
of this staircase generalises to 4d — 3 and 4d — 5. 

While in Sec. HIE this mapping will be employed to 
calculate diffusion constants during diffusive periods, we 
use it here to illustrate the concept of equivalence classes 
again. To this end, the mapping is generalised to incorpo- 
rate substrates as illustrated in shaded colours in Fig. 7: 
each substrate can be drawn as a box. This is seen as 
follows: because each spider being at a specific substrate 
site may either be in a spread or a relaxed configuration, 
a substrate has to be indicated at two different locations 
in the centre of mass space (thus the width of every box 
equals 2). Furthermore, since spider 1 being or not being 
at a substrate does not affect spider 2, every box indicat- 
ing a substrate at spider I's track has to be of a size that 
it contains all allowed configurations of spider 2, and vice 
versa. 

We now return to an example discussed in Sec. IIIB: 
Eq. (15) shows all configurations in which spider 1 has 
cleaved two more substrates than spider 2 and only spider 
2 is attached to a substrate. We referred to this set of 
configurations as the equivalence class [Acr]= [2i]. This 
situation is illustrated in Fig. 7, where there are A = 2 
more boxes (i.e. substrates) for spider 2 than for spider 1. 
The eight allowed configurations contained by the ellipse 
in this figure are only contained in one box (tr = 1), such 
that these states provide an geometrical interpretation of 
the equivalence class [2i] . 

Leaving the boundary period in this picture corre- 
sponds to removing the encircled box (i.e. cleaving the 
substrate) and stepping down (i.e. away from the sub- 
strate boundary). 



E. Diffusive Periods 

We now employ the mapping of the spider team mo- 
tion to a confined random walk in order to analyse the 
spider team's dynamics during a diffusive period. Let us 
first examine the transition rates between neighbouring 
points in the confined random walk picture. Consider, 
for example, the point 

(ci,c2) = (o,o) = ;°; (19) 

in Fig. 6. From this configuration, every leg may un- 
bind from its product with rate 1, and then rebind to ei- 
ther the same product site again, or move to the allowed 
neighbouring site at equal probability ^. In the confined 
random walk picture, this leads to transition rates of 1 • | 
along each connection between adjacent sites from (0, 0). 
The same argument applies to any site within the allowed 
region, so that the transition rate between any two lattice 
sites equals ^, cf. Fig. 8(a). This leads to the following 
master equation for the occupation probability Pci.c2 on 
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FIG. 7. (Colour online) Substrate in the staircase random 
walker picture (n = d = 2 as before). Like in Fig. 6, explicit 
configurations are shown for some points. In addition, boxes 
are drawn which correspond to the substrates on spider I's 
(vertical blue box), or spider 2's (horizontal red boxes) lane. 
This can be understood as follows: when a spider is attached 
to a substrate with its right leg, it can be either in the spread 
or the relaxed configuration. Hence a substrate at position c 
has to be indicated at two points in the centre of mass space, 
namely at c — | and c — 1; therefore the substrate boxes 
have width 2. Encircled in the figure are the eight states 
which have spider 2 at •oioo or o»i66, respectively, and 
spider 1 in one of the five states ••coo, oo»»6. The 
resulting states correspond clearly to those of Eq. (15) and 
Fig. 4. In the figure, there are three horizontal red boxes 
(substrates on spider 2's lane), and only one vertical blue box 
(substrate on lane 1). Hence, the difference of the product 
sea's ends is A = 2. Since the encircled states ® - ® have, 
by direct reading, only spider 2 at a substrate {i.e. they are 
only contained in cr = 1 box) , they form the equivalence class 
[A.l = f2il. 



the confined lattice: 

dt ''''' - 



(ci,C2> 



(20) 



where the sum runs over all nearest neighbours (ci,C2) 
of (ci,C2). In order to calculate the diffusion coefficient 
x'^{t)) we determine the time derivative 
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of the mean square displacement of the spider team: 



d 
dt 



(ci,c2)ec 



,2 V - 

'Cl,C2 / , 2 
(ci,C2) 



(ci,C2) 



(21) 

where Xc^^c2 — ^(ci + C2) is the position of the spider 
team on the molecular track for given values of ci and 
C2, and the summation extends over all {01,02) within 
the allowed region C. This equation can be reorganised 
such that 



dt 



[x\t)) 
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(22) 
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FIG. 8. (Colour online) Diffusion in the staircase environ- 
ment, (a) Transition rates between the sites of the staircase 
environment. Along every arrow drawn, the rate is | leading 
to local detailed balance, (b) The staircase can be split into 
elementary cells, numbered with integers. 



To evaluate this expression we split the lattice into el- 
ementary cells as shown in Fig. 8(b), and use that for 
asymptotically large times, t — > 00, the probability den- 
sity P varies only little between neighbouring elementary 
cells. This follows from translational symmetry; every 
cell obeys the same master equation. The master equa- 
tion, Eq. (20), then implies a nearly uniform probability 
distribution within each elementary cell j [33] . Upon as- 
suming a constant value Pj within each unit cell, carrying 
out the sum over an arbitrary elementary cell j leads to 
a further simplification 



c, 



(ci,C2> 



"(C1,C2) 



2 



3 1 



(23) 



independent of j. Altogether, we obtain 
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t->-oo dr 
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16 



where in (*) we used that each elementary cell comprises 
8 points, and in (f ) we employed the normalisation con- 
dition for P. This procedure can be generalised for arbi- 
trary d. The formula for the diffusion constants for n — 2 
then reads 

This theoretical result agrees well with simulation data 
for the diffusion constant _D, as a function of the leash 
length d, see Fig. 9. 



IV. OPTIMISATION OF DIRECTED MOTION 

In the previous sections we mainly focussed on en- 
semble properties of spider teams. However, in experi- 
ments or applications one has to deal with single reali- 
sations of the stochastic process, i.e. single trajectories. 
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FIG. 9. (Colour online) Diffusion constants as a function of 
the leash length d for n — 2 and 3 spiders. The dashed line 
shows the theoretical result for n — 2 (Eq. (25)); solid lines 
are asymptotics for d — ^ oo. Our theoretical approximation 
is in good agreement with simulation data (points). 



cf. Fig. 2(f). Since it is desirable to achieve a molecular 
motor design that works reliably, one would like to min- 
imise the randomness of the trajectory, i.e. the motion's 
standard deviation 



(26) 



It is interesting to ask how the microscopic properties of 
the spider team (n, d) influence a: can we optimise the 
performance of a spider team? Is there an optimal choice 
of parameters, n and d, which reduces the randomness of 
a spider teams' motion to a minimum? 

The randomness is determined by the interplay be- 
tween the dynamics of the spider team during its dif- 
ferent episodes of motion, i.e. the boundary periods and 
the diffusive periods. For each episode we found a char- 
acteristic feature: during boundary periods the spider 
team motion is essentially ballistic which can be quan- 
tified in terms of the mean number of consecutive steps 
(S) (cf. Eq. (17)). In contrast, during a diffusive period 
the spider team performs a random walk with a diffusion 
constant D (cf. Eq. (25)). 

We have already learned in Sec. IIIC and Fig. 5 that 
there is an optimal choice of parameters for the number 
of consecutive directed steps during a boundary period 
(see Fig. 5). One could now naively conclude that the 
predictability of a spider team motion can as well be 
optimised with the same set of parameters. However, 
this argument would overlook the impact of the diffusive 
periods. Indeed, there are several effects which influence 
the randomness during these episodes: 

(i) In Sec. HIE we noted that the diffusion constant D 
grows with the leash length d (Eq. (25) and Fig. 9). 
Since D determines the mean square displacement 
during a diffusive period, increasing d would then 
also imply a greater randomness, a. 



(ii) Conversely, a higher diffusion coefficient during dif- 
fusive periods speeds up all dynamic processes. 
Thus, in a given time window, larger d make it more 
probable for a spider team to return to the bound- 
ary and start moving ballistically [34]. 

The combined effect of these two processes can be esti- 
mated by analysing a random walker with an absorbing 
boundary. In one dimension, one finds that (^^(i)) 



\/ Dt [35, 36]. Hence, (i) and (ii) together would lead to 
an increase of a with d. 

(iii) Consider the geometrical interpretation of the tran- 
sition from the boundary period to the diffusive pe- 
riod as given in Fig. 7. In this picture, entering a 
diffusive period corresponds to removing the low- 
ermost red box, and stepping to one of the three 
points on the very bottom. Right after this transi- 
tion, the average minimal distance (xq) of the spider 
team from the boundary is therefore given by 



(27) 



as can be inferred from counting the different transi- 
tion pathways. Hence, with increasing d, the spider 
team entering the diffusive period remains closer to 
the boundary, and is thereby more likely to reenter 
a boundary period quickly. 

(iv) In Sec. HID we have shown that with increasing 
leash length d the number of pathways in state space 
to reenter a boundary period also increases. Pic- 
torially, this can be inferred from the mapping of 
the spider team's motion to a random walker in a 
staircase environment: the longer the leash length 
d the larger is the "angle" under which a random 
walker sees the boundary of the staircase. Thus, 
when the random walker takes an arbitrary direc- 
tion the probability that it walks toward the bound- 
ary is increasing with d. 

Since their is no unique trend in the various effects 
discussed above (i)-(iv), it is difficult to conclude what 
would be the dominant effect of the diffusive period on 
the randomness. Therefore, we numerically determine 
the randomness of the spider team during diffusive peri- 
ods [37]; this quantity is depicted in Fig. 10. We observe 
that the mean squared distance from the boundary is 
smaller for larger d at all times. This implies that — 
considering only diffusive periods — increasing d leads 
to a reduction of the randomness. From this we can infer 
that the effects (iii) and (iv), which decrease the random- 
ness of the process with increasing d, overcompensate the 
effects (i) and (ii). 

Altogether we can now conclude the influence of the 
diffusive periods as follows: 
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FIG. 10. (Colour online) Randomness during a diffusive pe- 
riod. Shown is the mean squared minimal distance to the 
boundary of a random walker in the staircase environment 
(n = 2), Fig. 6. The walkers start randomly along every 
point which provides an entrance to the diffusive period (for 
d — 2, these are the three points at the very bottom of Fig. 7); 
they are absorbed when they reach the boundary (which is the 
second substrate box in Fig. 7; note that the lowermost box 
has been removed when the walker entered the diffusive pe- 
riod!). Obviously, the mean squared distance is the greater 
the smaller d is. Increasing d thus decreases the randomness. 

Analogously we can decipher the influence of boundary 
periods. Going back to Fig. 5 we observe: 

d<d°f^: d/' ^ (S) /■ ^ <7\, 
d>d°f^: d/" ^ {S)\ ^ a /- . 

These considerations explain that i/ there is an opti- 
mal value , it must be found beyond . This is in 
agreement with our data: Figure 11 shows the existence 
of a minimum of the randomness, and its positioning with 
respect to . Remarkably, the positions of both optima 
are strongly correlated (see Tab. II). 

In conclusion, our analysis shows that the randomness 
of the spider team is mainly determined by the mean 
number of steps (S") during boundary periods. Diffusive 
periods have only a small effect on the randomness and 
change the optimal parameters only slightly. 

V. CONCLUSIONS 

Based on existing models for molecular spiders [11, 12], 
we proposed a model for a spider team that explores the 
collective behaviour of cooperating spiders: in our model, 
bipedal spiders are jointly attached to a (zero-mass) link- 
ing cargo. Each spider walks on its own one-dimensional 
track. This leads to a spacial constraint which can be 
characterised by the maximal span d, of the resulting 
spider team. 



FIG. 11. (Golour online) Standard deviation a of the spiders 
movement and the mean number of steps {S) as a function 
of the leash length d. Both a and (5*) show extrema. To em- 
phasise the correspondence between the minimum of a and 
the maximum of (S) (cf. Fig. 5), the (S)-axis is drawn in 
reverse (see right scale), a is measured at the time t* when 
the mean displacement {x) equals 1000. This choice is arbi- 
trary; for smaller values the minima of a persist, but are less 
pronounced, cf. Fig. 2(b). 



Depending on the cleavage rate of the substrate r < 1, 
the number of coupled spiders n, and the leash length d, 
we found that the coupling leads to a significant enhance- 
ment of many of the spider's motor properties: spider 
teams show a significant increase of their mean displace- 
ment; their motion is a lot less random; the ensemble's 
velocity can be increased by more than an order of mag- 
nitude; and the superdiffusive behaviour lasts longer for 
orders of magnitude in time. Unlike single spiders, co- 
operating spiders could therefore — at least in theory — 
be employed for executing well-defined tasks reliably. 

Like their individual counterparts [16], spider teams' 
motion can be characterised as being in either a bound- 
ary, or a diffusive period. We found that the charac- 
teristic quantity is the mean number of consecutive di- 





n = 2 




n = 3 


r 


,opt 


dT 


jopt ,opt 


0.001 
0.01 
0.02 
0.05 
0.1 
0.2 


~10-11 
4 
3 
2 
2 
2 


~13 
5 
4 
3 
3 
2 


5 6 
3-4 4-5 

2 3 
2 3 


TABLE II. Optimal values of d for n = 2 and 3, and several 
values of r. Compared are the values of d which maximise 
the mean number of steps during a boundary period, and 



that which minimise the randomness (for a comparison see 
Fig. 11). Clearly, both values of d are closely correlated, where 
ci?g^ is only slightly smaller than d^^'. 
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rected steps, (S), which a spider team performs during a 
boundary period. In simphfied language, (S) integrates 
the number of steps which the spider teams walks, as 
long as it stays in the vicinity of the co-moving bound- 
ary between substrate and product sites. (S) is closely 
related to the bias p+ of single spiders [12]. For small 
r, we succeeded in calculating (S) analytically through 
an equivalence class formalism which made use of the 
time scale separation of dwell times on products and 
substrates. This formalism is exact for r — > 0, regard- 
less of the number of coupled spiders and the tightness 
of the coupling. We explicitly calculated values for vari- 
ous small spider teams, and find excellent agreement with 
simulation data. For small but finite cleavage rates r, the 
formalism still holds as an approximation for relatively 
tight coupling. We found that in this case there is an op- 
timal value for the coupling tightness d which maximises 
the mean number of steps. 

Next, we provided a mapping of the stochastic mo- 
tion of an n-spider team to a random walker in an n- 
dimensional environment. The motion is confined be- 
tween two parallel boundaries which have the shape of 
staircases. This mapping is exact and allows a comple- 
mentary interpretation for the equivalence classes: sub- 
strates can be drawn as boxes which are easy to enter 
for random walkers but impossible to leave without re- 
moving, which happens slowly on a timescale r~^. It is 
then straightforward to see that an equivalence class cor- 
responds to an intersection of boxes (cf. Fig. 6). The 
staircase picture also allows to quantify the dynamics 
during the diffusive periods of spider teams: in that case, 
boxes can be ignored and spider teams correspond to or- 
dinary diffusive random walkers on the confined lattice. 
We calculated the diffusion constants for 2-spider teams 
and find good agreement with simulation data. 

The analysis of the mean number of consecutive steps 
during a boundary period, (S) (which shows a maximum 
for some value of the leash length d) , taken together with 
the diffusion constants D (which grow with d) allow for a 
comprehensive explanation of our observations. We show 
that the optimal value of d that minimises the random- 
ness (which involves boundary and diffusive periods) dif- 
fers only slightly from the leash length maximising the 
mean number of steps during a boundary period; see 
Fig. 11. 

The staircase picture also illustrates that despite the 
difference in complexity, a single spider and a spider 
team can both be described by similar effective random 
walk models: the motion of a bipedal spider which has 
a non-trivial stepping gait can be fully described by its 
centre of mass coordinate which performs simple one- 
dimensional random walks [12]. Likewise, the motion of 
an n-spider team which involves complicated interactions 
between the spiders can equivalently be described by 
another single coordinate which performs n-dimensional 
random walks that are however geometrically confined 
due to the leash constraint. 

Our results show that the primary factor for improving 



the motor properties of molecular spiders is the accessi- 
bility of substrate sites for the spider legs: while single 
spiders only have access to one substrate at a time, an 
n-spider team can reach n substrates. This would im- 
ply that there is a significant difference between truly 
one-dimensional spiders [12] and quasi one-dimensional 
spiders [10]. This is enforced by a very recent study of 
Olah et al. [19] who examined molecular spiders on a nar- 
row 2-dimensional lattice. As well, it is in full accordance 
with recent data by Samii et al. [18] who concentrated on 
hand- over- hand spiders: they showed that motor proper- 
ties of this class of spiders which have access to more than 
one substrate site at a time are superior to inchworm spi- 
ders which can only reach one substrate at once [17, 18]. 

The results presented here can be extended in multi- 
ple ways. In analogy to individual spiders, further stud- 
ies could concentrate on varying design specifics like the 
number or the length of legs [18]. Likewise, the un- 
derlying chemical processes [9, 10] could be modelled in 
greater molecular detail also for spider teams. Similarly, 
the team's spiders' stepping gait could be varied, poten- 
tially profiting from studies about the motion of individ- 
ual hand-over-hand spiders with more than two legs [18] 
which seem to be difficult to realise in the experiment. 

Unlike other studies (e.g. [16]) which have extensively 
investigated the role of the cleavage rate r, our focus was 
different and the variation of r was only a side aspect 
of this work. Nevertheless, our analysis hints towards a 
scaling behaviour which maps the quantity {S{r,d)) to a 
universal form S{d) which is independent of r. In this 
spirit, it would also be interesting to study the connec- 
tion of the optimal leash length and the cleavage rate r. 
It appears that this relation might be rather simple for a 
wide parameter range, although its mathematical formu- 
lation seems to be very complex. The difficulty is that 
the simplified formulation of the problem presented here, 
i.e. the equivalence classes, can not be applied directly. 
One possibility to address this problem might lie in draw- 
ing analogies from related models like the burnt-bridge 
model [38]. For example, it has been studied for dimeric 
motor molecules [39] and as an exclusion process [40]. 

Our results might also be relevant to study collective 
properties of molecular motor assemblies theoretically, cf. 
Ref. [41] and the references therein. These models are 
relevant to understand the interplay between biological 
motor molecules like kinesin, dynein and myosin inside 
cells [42, 43]. In contrast to spiders, biological motors 
are fuelled by ATP hydrolysis; they can build up signif- 
icant pulling forces due to strong mechanochemical cou- 
pling [44]. In particular, recent experiments addressed 
the complex interplay of multiple coupled kinesin mo- 
tor proteins where the motors are coupled via a DNA 
leash of certain length. It is interesting to note how in 
these experiments teams of two kinesin motors outper- 
form a single motor in terms of run-length and pulling 
forces [45-47] . Similarly, cooperative effects also improve 
the properties of two coupled burnt-bridge motors mod- 
eling coUagenase transport [39]. 
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In conclusion, we believe that our model of coupling 
molecular spiders provides insight on how cooperative be- 
haviour evolves on the molecular scale. We hope that our 
ideas about molecular spiders help advance a young and 
fast growing field in which much focus is put on the con- 
struction of novel, more efficient, molecular designs [4]. 
We believe that our findings are not limited to the case 
of molecular spiders, but apply to molecular machines 
working together in general. 



fusive period. These read 

P(J|[0])=MJ-2I[1]), 

p(jIW)-^(pO--^|[*-i])+p(j-^|[* + i])) 

pij\[d-l])^\[pij-\\[d-2])+pij\[d])), 

pij\[d])=Up{j^l\[d^l]), 

where < i < d — 1. The mean number of steps {S{x)) 
which a spider team walks from class [x] until going to 
the diffusive period is then given by 
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Appendix A: Derivation of Eq. (17) 



We analyse the graph for a 2-spider team with arbi- 
trary d as depicted in Eq. (14). According to this graph, 
transitions [z] — !■ [i± l] are equally likely as long as i < d, 
whereas [d] — >■ [d — l] happens at probability 11. During 
every transition, the spider team performs a fractional 
step \ = \- Only during the transition [d] — > [d — l] , 
no step is integrated; in return, [d — l] — >■ [d] leads to 
a whole step for the team. This is due to the very defi- 
nition of the number of steps during a boundary period, 
which comprises all cleavages but for each spider's last 
cleavage before the team enters the diffusive period. 

With these preparations, we can now establish the 
probabilities p(i|[i]) that a spider team, being in class 
[z] , performs exactly j steps before leaving into the dif- 



{S{x))= ^p^^\^' 

J=0,i,... 



(A2) 



Inserting Eq. (A2) into Eq. (Al), and by renumbering 
indexes we obtain 

{sm = \ + {s{i)) . 

(^W> = ^ + 2(^(^-1)) + + (A3) 
(5(d-l)) = i + i(5(d-2)) + i(5(d)), 

(5(d)) = n + n(5(d-i)), 

where again < i < d — 1. Solving this system of equa- 
tions, we obtain the recursion relation 



(5(fc)) = (5(fc + l)) + fc 



1 



(A4) 



for < fc < d — 2. Substituting this into the remaining 
equations leads to 



(A5) 



and finally 



(5(d-i))^^-iH^-i)r^ + T^.(A6) 

Since a spider always enters a boundary period in class 
[d— 1] in the limit r ^ (cf. Eq. (13)), the last equation 
is equivalent to (S*), Eq. (17). 
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